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Estimates  of  the  rotational  distributions  in  the  N2(C  n^,  v'  = 

0,1)  product  of  the  reaction 

Ar*(3p2^Q)  +  N2(X  ^Eg)  — Ar(ls)  +  N*(C  \) 

have  been  made  on  the  band  profiles  of  the  HJC  \^,   v'  =  0.1)  -^ 
N2(B  Ilg,  v"  =  0)  fluorescence.  These  measurements  were  made  at  relative 
energies  of  0.074  eV,  0.089  eV,  and  0.161  eV.  The  N2(C  ^n^^,  v'  =  0)  band 
was  found  to  be  well  represented  by  Boltzmann  distributions  with  charac- 
teristic temperatures  of  1700  +  100  K,  1600  +  100  K,  and  2200  +  100  K, 
respectively.  The  N2(C  "^n^^,  v'  =  1)  data  could  not  be  fit  by  Boltzmann 
distributions,  but  did  fit  "Golden  Rule"  distributions  within  the  experi- 
mental error.  The  reaction  involving  metastable  argon  and  ground  state 
hydrogen  molecules  was  examined  by  observing  a  continuum  emission.  The 
cross  section  for  fluorescence  was  measured  and  seen  to  rise  rapidly  from 
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an  onset  at  0.080  eV  and  reach  a  maximum  around  0.125  eV.  The  spectral 
distributions  were  seen  to  peak  at  lower  energies  than  that  calculated 
by  James  and  Coolidge,  leading  to  speculation  that  at  least  some  of 
the  photon  emission  may  arise  from  an  excited  state  of  ArH  rather  than 
the  excited  Ho. 
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CHAPTER  I 
INTRODUCTION 


A.  Molecular  Beams 

The  first  detection  of  a  molecular  beam  was  Fleming's  discovery 
of  a  shadow  on  the  glass  walls  of  incandescent  lamps  with  copper  fila- 
ments [1].  At  the  time,  this  was  considered  a  confirmation  of  the  kinetic 
theory  which  predicted  that  a  molecule  would  travel  in  a  straight  path 
in  the  absence  of  collisions  or  applied  fields.  This  phenomenon  was 
previously  observed  for  ions  and  electrons,  but  this  was  the  first 
discovery  involving  neutral  species. 

The  first  deliberately  designed  experiments  utilizing  molecular 
beams  were  the  experiments  of  Dunoyer  [2],  which  again  verified  the  idea 
that  a  molecule  which  is  not  disturbed  by  collisions  will  continue  to 
travel  in  a  straight  line. 

The  next  significant  work  involving  the  use  of  molecular  beams 
was  the  work  of  Stern  which  ultimately  led  to  the  famous  Stern-Gerlach 
experiments  verifying  the  existence  of  space  quantization  and  electron 
spin  [3].  This  work  led  to  techniques  for  the  measurement  of  the  mag- 
netic moment  of  various  molecules  by  the  magnetic  resonance  techniques 
of  Rabi  [4],  and  the  further  refinements  in  molecular  beam  techniques 
by  Ramsey  [5]. 

During  this  time,  the  study  of  chemical  reactions  was  moving  from 
the  field  of  bulk  reaction  studies  toward  attempts  to  understand  the 

1 


dynamics  of  such  reactions  in  terms  of  individual  collisions  between 
molecules.  It  has  been  in  this  field  that  the  molecular  beam  has  found 
its  most  significant  uses.  The  pioneering  experiments  in  this  field  were 
conducted  with  effusive  or  oven  beams,  which  provided  only  thermal 
averages  of  the  quantities  being  measured,  or  else  were  velocity  selected, 
greatly  reducing  the  beam  intensity.  The  first  major  improvement  in 
this  condition  was  the  development  in  the  1950s  of  the  multiparallel  canal 
sources.  In  1953,  King  and  Zacharias  [6]  used  a  bundle  of  hypodermic 
needle  tubing  and  sheets  of  foil  to  produce  a  beam  of  high  intensity. 
This  was  followed  rapidly,  in  1955,  by  the  production  of  intense  beams 
of  anmonia  by  the  multiparallel  canal  source  for  the  study  of  masers  17 J. 

At  the  same  time  the  first  nozzle  beam  sources  were  being  devel- 
oped for  use  in  studies  of  chemical  interest.  In  1954  Becker  and  Bier 
[8]  first  achieved  the  high  pumping  rates  required  for  a  successful  nozzle 
beam.  Because  of  the  higher  intensities  and  narrower  velocity  distri- 
butions of  this  type  of  beam,  the  nozzle  beam  has  shown  a  marked  increase 
in  use  as  high-speed  vacuum  systems  have  come  into  more  common  use. 

B.  Supersonic  Nozzle  Beams 

The  supersonic  nozzle  beam  is  a  source  of  highly  collimated,  high- 
intensity,  and  nearly  monoenergetic  particles.  The  beam  is  created  by 
applying  a  high  pressure  behind  a  small  orifice  and  allowing  the  gas  to 
expand  into  a  vacuum.  The  expansion  jet  is  shown  in  Figure  I-l. 

The  major  feature  of  interest  is  the  central  region  terminated  by 
the  Mach  disk.  Within  this  region  is  the  supersonic  portion  of  the  expan- 
sion. By  the  interception  of  this  portion  of  the  gas,  a  supersonic  beam 
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produced.  This  interception  is  accomplished  by  means  of  an  appropriately 
positioned  device  called  a  skimmer.  The  skimmer  penetrates  into  the 
Mach  disk,  and  shields  the  supersonic  flow  from  collisions  with  the 
background  gas.  These  collisions  would,  if  not  prevented,  destroy  the 
desired  beam  characteristics.  Beyond  the  skimmer  is  a  low-pressure 
region  in  which  the  beam  may  propagate  without  this  interference. 

As  the  beam  is  essentially  formed  at  the  mouth  of  the  skimmer,  the 
properties  of  the  gas  at  this  point,  the  gas  pressure  and  velocity  distri- 
bution, determine  the  ultimate  characteristics  of  the  beam.  These  condi- 
tions are  described  by  [9] 

dn  =  n(m/2  Kt)^/^exp  -(m/2kT)(c-u)^  d^c  (1) 

where  n,  T,  and  u  are  the  density,  translational  temperature,  and  mean 
gas  velocity.  The  mass  of  the  molecule  in  the  beam  (m)  and  the  Boltzmann 
constant  (k)  are  the  remaining  parameters  in  this  equation.  From  inte- 
gration of  Equation  (1),  several  important  relationships  are  derived. 
These  are  given  as  Equations  (2)  through  (7). 

X^/6*   =  0.67(Pq/P^)^/^  (2) 

u  =  ((2/(y  -  l))(kT^/m)^^2  ^3^ 

M  =  u/(YkT/m)^/2  (4) 

S  =  u/(2kT/m)^/2  =  (y/2)^^^M  (5) 

1(6=0)  =  (ATniao/27r)l^)M(3  +  yM^)[1  +  (y  -  ^)^'^/2]''^^^       (6) 


In  the  above  equations  the  symbols  are  M  =  Mach  number  =  the  ratio  of 
the  flow  velocity  to  the  local  speed  of  sound,  u  =  the  bulk  stream  veloc- 
ity, k  =  the  Boltzmann  constant,  d*  =  orifice  diameter  of  the  nozzle, 
Y  =  specific  heat  ratio  (C  /C  ),  T  =  temperature,  m  =  mass,  p  =  pressure, 
S  =  speed  ratio  =  the  ratio  of  velocity  of  the  beam  to  that  of  a  thermal 
beam,  n  =  number  density,  A  =  area,  I  =   distance  downstream  from  the 

skimmer,  G  =  mass  flow  through  the  nozzle  in  grams/second,  a  =  speed 

1/2 
of  sound  =(YkT/m)  '  ,  x  =  distance.  The  subscripts  are  0  =  stagnation 

(source)  conditions,  1  =  nozzle  exhaust  chamber,  skimmer  orifice,  2  = 
collimation  chamber,  collimator  slit.  The  above  equations  are  quite  use- 
ful in  design  considerations  for  such  a  source. 

Another  key  feature  of  the  supersonic  beam  is  a  "cooling"  of  the 
internal  degrees  of  freedom  providing  additional  translational  energy  to 
the  beam.  This  follows  from  the  thermodynamics  of  an  ideal  gas  under- 
going isentropic  expansion.  The  total  amount  of  energy  available  to  the 
gas  is  given  by 

E  =  CpTg  (8) 

where  C  is  the  heat  capacity  at  constant  pressure  of  the  gas  and  Tq  is 
the  temperature  of  the  gas  behind  the  nozzle  (stagnation  temperature). 
Recalling  that  C  =  y/(y-1)»  (y  =  C  /C  ).  If  not  all  of  this  energy  is 
converted  to  translational  energy  the  beam  is  characterized  by  a 
"temperature"  T  and  the  energy  converted  is  given  by 

E  (Y/(Y-l))(Tn-T)  =  mu^/2  (9) 


where  u  is  the  velocity  of  the  beam.  The  ratio  of  the  stagnation 
temperature  to  the  final  beam  temperature  is  given  as 

Tq/T  =  1  +  ((y  -  1)/y)M^  (10) 
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where  the  Mach  number  (M)  is  defined  as  u/lykT/M)  '  . 

This  temperature  may  be  broken  down  into  components  relating  to 
random  translational  motions,  vibrational  energies,  and  rotational  ener- 
gies. By  assuming  most  of  the  rotational  energy  relaxation  occurs  in  the 
high-pressure  portion  of  the  expension  where  collisions  between  beam 
members  predominate,  it  is  possible  to  calculate  a  rotational  temperature 
for  a  given  expansion  [10].  This  assumes  the  collision  process  is 
described  by  a  van  der  Waals  type  of  interaction  potential  (V  =  -C5/r^). 
The  results  are  given  as 

where  t(y,^)  is  a  tabulated  constant,  a  measure  of  the  strength  of  the 
rotational  to  translational  energy  coupling,  y  is  the  ratio  C  /C^, 


Xq  =  (3/2)2/3(3  Cg/kT^)^/^  (5/3)  (n  /2)(32/27)  sin^  9  (3/2)M 

(12) 


and 


Zq  =  r2  K[{(y  -  1)/(Y  +  l)}^/2  (2/Y  +1)^/^"'']  (13) 


where  k  is  the  BoUzmann  constant,  n^  is  the  number  density  of  the  gas 
at  the  nozzle,  the  geometric  term,  sin^  e,  is  calculated  to  be  0,0812, 
Tq  is  the  nozzle  aperture  radius,  and  K  is  a  peaking  factor  related  to  y. 

C.  Capillary  Array  Beams 

Another  type  of  molecular  beam  source  which  has  become  popular 
is  the  multichannel  or  capillary  array  source.  Although  its  centerline 
intensity  is  not  as  great  as  a  nozzle  beam  source,  the  capillary  array 
beam  is  a  high- intensity  directional  molecular  beam  source  when  compared 
with  the  conventional  oven  source.  The  capillary  array  source  has  the 
advantages  over  a  nozzle  source  of  a  much  lower  gas  flow  rate  and  operat- 
ing conditions  which  may  be  accurately  recreated,  allowing  for  greater 
day-to-day  reproducibility  of  the  beam. 

Pauly  and  Toennies  [11]  have  calculated  a  theoretical  centerline 
intensity  for  a  conventional  oven  source  operating  under  Knudsen  flow 
conditions,  that  is,  the  mean  free  path  of  the  particle  in  the  source 
is  on  the  order  of  the  exit  hole  diameter.  This  relationship  is 

1(9  =  0)  =  1.12  X  10^^  PqFq/(MT)^/^  (14) 

where  Pq  is  the  source  pressure  in  torr,  Fq  is  the  orifice  area  in  cm  , 
T  is  the  temperature  in  degrees  Kelvin,  and  M  is  the  molecular  weight 
of  the  gas. 

When  this  is  compared  to  a  multichannel  array  in  which  the  length 
of  the  channel  is  much  larger  than  the  diameter  of  the  channel  and  which 
is  operated  under  conditions  such  that  the  mean  free  path  is  on  the  order 
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of  the  channel  length,  a  relationship  between  the  centerline  intensities 
of  the  two  sources  may  be  obtained  as 

^(^  =  °Wray/^^Q  =  °)oven  =  ^'^^  V^/^CFtm)^/^  /(J/^)   (15) 

where  V  is  the  average  velocity  in  the  source,  N  is  the  total  gas  flow 
rate,  a  is  the  scattering  cross  section,  F  is  the  total  area  of  the 
source,  t  is  the  transparency  of  the  array,  and  m  is  the  total  number  of 
holes  in  the  array.  Becker  and  Houkes  [12]  have  demonstrated  that  for 
a  given  beam  intensity  the  gas  flow  from  a  capillary  array  source  is 
approximately  100  times  less  than  the  flow  from  a  conventional  oven 
source,  thereby  reducing  the  stress  placed  on  the  vacuum-pumping  system 
and  conserving  what  may  be  an  expensive  quantity  of  gas. 

D.  Reaction  Cross  Sections 
At  a  particular  reaction  energy,  an  expression  for  the  reaction 
rate  (R),  as  determined  for  crossed  molecular  beam  experiments,  may  be 
written  as 

R  =  a(v^)v///n^(x,y,z)ng(x,y,z)  ^^  (16) 

where  a(v  )  is  the  reaction  cross  section,  v  is  the  relative  velocity  of 
the  collision  partners,  n»(x,y,z)  and  ng(x,y,z)  are  the  coordinate 
dependent  beam  densities  of  the  two  beams,  and  V  is  the  interaction  vol- 
ume. The  primary  experimental  difficulty  in  using  this  equation  is  the 
measurement  of  the  quantities  n.  and  ng.  Assuming  these  are  at  most 
slowly  varying  quantities  within  the  interaction  region,  the  reaction 
rate  may  be  written  as 


where  the  beam  densities  are  no  longer  spatially  dependent.  In  experi- 
ments such  as  the  one  to  be  described,  one  of  the  direct  products  of  a 
collision  is  a  photon  which  is  later  detected.  This  photon  flux  is  then 
proportional  to  the  reaction  rate. 

Sp  =  KR  =  Ka(v^)v^n^ng  (18) 

where  S  is  the  photon  flux  and  K  is  the  proportionality  constant.  The 
metastable  beam  detector  measures  the  intensity  (I^  =  v^/n^^)-  From  this 
an  expression  for  the  reaction  cross  section  may  be  obtained: 

a(v^)  =  K-^SpV^/VgV^  (19) 

Upon  integration  over  all  collision  energies  with  a  Maxwell -Boltzmann 
weighting  factor  this  quantity  may  be  related  to  the  more  familiar  rate 
constant,  which  is  determined  in  most  kinetics  experiments. 

E.  The  Franck-Condon  Principle 
In  1925  Franck  [13]  proposed  a  method  for  explaining  the  band 
spectrum  observed  in  diatomic  molecules.  His  basic  postulate  was  that 
the  electron  transition  directly  affects  neither  the  position  nor  the 
momentum  of  the  nuclei.  That  is,  during  an  electronic  transition,  the 
electron  will  jump  from  one  potential  surface  to  another  more  rapidly 
than  the  nuclei  can  respond  to  the  impetus  imparted. 

In  1927  Condon  [14]  made  the  first  attempts  to  connect  this  pos- 
tulate with  quantum  mechanics.  In  1928  Condon  [15]  made  a  more  elaborate 
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proposal  on  this  relationship.  Using  the  notation  of  this  paper,  the 
mathematical  formulation  takes  the  following  form. 

Beginning  with  the  Born-Oppenheimer  approximation,  it  is  possible 
to  write  the  energy  of  a  state  E(e,n)  characterized  by  e,  an  ensemble  of 
all  the  electronic  quantum  numbers,  and  by  n,  the  vibrational  quantum 
number  (for  the  purposes  of  this  study,  rotational  motions  are  neglected), 

E(e,n)  =  Eg  +  E^  (20) 

Also  the  total  wave  function  may  be  written  as 

^en^r,x)   =  il)Q{X)rl)j^{r)  (21) 

Any  given  transition  (e',  n')  -»-(e",n")  may  now  be  given  in  terms 
of  i(;g^(x,r)  and  the  electronic  moment  M(e' ,e",n' ,n").  This  transition 
moment  may  be  resolved  into  components  depending  only  on  the  individual 
coordinates  as  Mg  +  M  .  The  transition  moment  R  may  now  be  calculated  as 

R  =  /Meil^^XW  dx  + /M^il^^X'i'eK  ^^  (22) 

As  M  does  not  depend  on  the  electronic  coordinates,  the  second  integral 
may  be  written  as 

/M^K^X  d-z^rr^*^'^   dTg  (23) 

Because  of  orthogonality  of  the  electronic  states  the  integral  reduces 
to  zero.  Thus  the  transition  moment  reduces  to 

R  =  /4^;ii>;;  dT^m^Ve*^'^  ^Tg  (24) 


n 


The  latter  portion  of  this  integral  is  the  electronic  transition  moment 
(R  )  and  in  the  Franck-Condon  limit  is  assumed  to  vary  only  yery   slowly 
with  nuclear  distance  and  as  such  may  be  replaced  by  an  average  value 
of  the  electronic  transition  moment  (R  ).  The  transition  moment  may 
then  be  written 


R  =  Re  /'f'^;  dr 


and  the  intensity  of  a  given  transition  may  be  written 


I    =  (64/3)77  cNy.v  Rg 


fKK  ^"^ 


(26) 


The  key  result  of  this  manipulation  is  that  the  intensity  of  the 
transition  is  proportional  to  the  square  of  the  overlap  integral  between 
the  vibrational  states  involved.  Therefore,  the  maximum  will  be  found 
for  a  transition  for  which  this  factor  is  maximized,  rather  than  for  a 
straightforward  vertical  transition. 

F.  Energy  Distributions  in  Product  States 
The  most  familiar  distribution  of  energy  in  a  product  is  observed 
in  bulk  reactions.  The  energy  in  this  process  has  had  time  to  randomize 
into  the  well-known  thermal  or  Boltzmann  distribution.  This  distribution 
is  characterized  by  a  well-defined  temperature. 

A  second  way  of  dividing  the  energy  in  the  product  states  is  given 
by  the  Fermi  Golden  Rule  [16].  In  this  case,  the  relative  transition 
probabilities  for  the  formation  of  a  given  final  state  is  given  by 
first-order  perturbation  theory  as 
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W^.^f  =  (27r/1i)|<f|V|i>|^  p^(e)  (27) 


where  V  is  the  interaction  potential  connecting  the  initial  (li>)  and 
final  (<f|)  states,  and  p^(e)  is  the  density  of  states  term  which  is 
dependent  on  the  available  energy  (e)  [17].  In  the  limit  of  a  highly 
impulsive  or  sudden  collision,  the  matrix  elements  in  Equation  27  become 
the  Franck-Condon  factors  connecting  the  initial  and  final  states 

|<f|V|i>|^-|<f|i>l2  (28) 

For  a  vibrating  rotor,  the  density  of  states  function  is  given  by 
(18) 

Pf(e)-(1  -  fv.)^^yi(l  -  %>)^^^  (29) 

where  f  ,  is  the  fraction  of  the  mean  available  energy  channeled  into 
a  particular  vibronic  state. 

Kinsey  [19]  has  calculated  a  density  of  states  function  for  a 
variety  of  conditions  including  the  partial  resolution  of  product  and 
rotational  states.  For  the  case  in  which  the  translational  energy  (E  ), 
total  energy  (E),  vibrational  state  (v),  and  rotational  state  (J)  are^ 
known,  the  density  function  is  given  by 

p^,(E',v',J')  =  (2J'  +  l)p^[E'  -E'(v',J')]  (30) 

where  the  density  of  translational  states  (p  )  is  given  by 

P,(E^)  =  .2"^"h]'^-^  -   A^Ey^  (31) 
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where  y  is  the  reduced  mass  of  the  colliding  molecules.  For  the  purpose 

of  calculating  the  internal  energy  of  the  product  states  the  rigid-rotor- 

harmonic-oscillator  model  has  been  used  with  first-order  correction  for 
anharmonicity.  This  gives  an  internal  energy  of 

Ej(v,J)  =  (v  +  l/2)(x)g  +  (v  +  l/2)^WgXe  "^  ^e^^^-^^) 

+  ag(v+  1/2)J(J+1)  (32) 

G.  Electronic  Energy  Transfer  Processes 
Electronic  energy  transfer  processes  are  observed  both  in  nature 
and  in  man-made  laboratory  experiments.  In  the  rarefied  gases  of  the 
upper  atmosphere,  electronically  excited  species  play  an  important  role 
in  the  chemistry  of  this  region.  Because  of  the  low  pressures,  and  cor- 
responding long  path  lengths,  energy  carriers  must  be  very  efficient  at 
transferring  energy;  that  is,  the  transfer  must  be  accomplished  in  a  wery 
few  collisions.  It  has  been  shown  that  vibrational  energy  often  requires 
tens  of  thousands  of  collisions  to  transfer  energy  [20],  while  energy 
stored  in  electronic  energy  levels  is  transferred  after  only  a  very  few 
collisions  [21]. 

In  the  laboratory  one  of  the  more  notable  successes  of  the  study 
of  electronic  excitation  processes  is  found  in  various  laser  systems, 
such  as  the  He-Ne  laser  [22],  CO2  laser  [23],  and  the  various  eximer 
lasers  [24].  Although  these  studies  have  led  to  many  practical  benefits, 
the  study  of  electronic  energy  transfer  has  also  led  to  greater  under- 
standing of  the  interactions  of  atoms  and  molecules. 
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Electronic  energy  transfer  takes  place  in  a  variety  of  processes, 
depending  upon  the  energy .available  and  the  dynamics  of  the  collisions. 
Some  of  the  more  common  processes  are  listed  in  Equations  (33)  through 
(37)  below. 

A*  +  XY  -.»  A  +  XY"^  +  e"  (33) 

(Penning  Ionization) 

A*  +  XY  -^  AXY"^  +  e"  (34) 

(Associative  Ionization) 

A*  +  XY  -*>  A  +  X  +  Y"^  +  e'  (35) 

(Dissociative  Ionization) 

A*  +  XY  -^  A  +  X  +  Y*  (36) 

A*  +  XY  -^  A  +  XY*  (37) 

In  all  of  these  reactions,  the  asterisk  indicates  the  location  of  the 
electronic  excitation.  In  the  last  two  equations,  (35)  and  (36),  the 
excited  product  will  relax  to  lower  states  by  the  emission  of  a  photon  if 
the  optical  transition  is  allowed. 

The  most  prolific  experimental  technique  in  the  study  of  neutral 
products  has  been  the  flowing  afterglow  technique  used  by  Setser,  Stedman 
and  Coxon  [25].  This  method  has  been  used  to  study  many  reaction  systems, 
obtaining  total  cross  sections  as  well  as  relative  cross  sections  for 
energy  transfer  into  various  rotational,  vibrational,  and  electronic 
states.  However,  this  type  of  experiment  yields  only  thermal  averages  of 
the  measured  quantities.  For  the  energy  dependence  of  the  processes. 
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molecular  beams  have  played  an  important  role.  The  technique  of  molecu- 
lar beams  has  been  used  to  study  this  type  of  process  in  time-of-flight 
crossed  beam  measurements  [26],  angular  scattering  measurements  [27], 
and  crossed  beams  in  which  photon  emission  is  used  as  the  detection 
process  [28].  It  is  the  latter  of  these  techniques  which  is  employed  in 
this  study. 

H.  Purpose  and  Scope  of  Present  Study 
The  reaction 


Ar*(%^Q)  +  n^ixh*)  -^  Ar(^S)  +  H^{c\) 


has  been  studied  by  the  detection  of  fluorescence  from  the  process 


N2(c\)  ^  N2(b\  )  +  hv 


by  Sanders,  Schweid,  Weiss  and  Muschlitz  [29]  who  determined  the  cross- 
section  response  to  collision  energy  for  the  total  radiation  emitted. 
This  process  has  also  been  investigated  by  Lee  and  Martin  [30],  using 
time-of-flight,  with  differing  results  for  the  onset  energy  for  the 
process.  Cutshall  and  Muschlitz  [31]  have  studied  the  energy  dependence 
of  the  distribution  of  the  energy  into  the  various  vibrational  levels  of 
the  N2{C)  manifold.  The  results  of  those  experiments  have  been  explained 
theoretically  by  Gislason,  Kleyn  and  Lus  [32]  using  a  model  in  which,  as 
the  reactants  approach  a  critical  distance,  the  NpCX)  potential  surface  is 
disturbed  by  the  close  lying  N^  +  Ar~  potential  energy  surface.  This 
causes  the  N^  molecule  to  begin  to  vibrate.  Several  vibrational  periods 
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later,  the  collision  partners  reach  a  second  critical  distance  of  approach 
and  jump  from  the  N2(X)  +  Ar*  surface  to  the  N|(C)  +  Ar(  S)  potential  sur- 
face. It  is  this  time  interval  prior  to  the  final  surface  jump  that  then 
influences  the  final  distribution  of  vibrational  states.  The  first  por- 
tion of  this  work  consists  of  an  investigation  of  the  rotational  distri- 
bution of  the  energy  within  specific  vibrational  levels. 

Investigations  of  the  radiation  from  the  hydrogen  continuum  have 
been  carried  out  theoretically  [33,34,35,36]  and  experimentally  by 
Coolidge  [37]  as  well  as  by  Smith  [38]  and  by  Finkelnburg  and  Weitzel 
[39].  The  potential  energy  for  the  Hp  molecule  has  been  recalculated 
recently  by  Kolos  and  Wolniewicz  [40].  The  second  portion  of  this  work 
consists  of  measurement  of  the  Ar*/H2  cross  section  for  the  observed 
continuum  radiation. 

Chapter  II  is  a  description  of  the  apparatus,  the  vacuum  system, 
and  the  data  collection  system.  Chapter  III  presents  a  description  of 
the  general  experimental  procedures  as  well  as  details  of  the  specific 
systems  studied.  Chapter  IV  presents  the  measured  spectra  and  cross 
sections  as  well  as  a  discussion  of  the  errors  in  the  data.  Chapter  V 
is  a  discussion  of  the  results  of  this  study.  Previous  measurements 
are  discussed  when  available,  and  possible  mechanisms  for  these  processes 
are  presented. 


CHAPTER  II 
DESCRIPTION  OF  THE  APPARATUS 


A.  Introduction 

The  experimental  apparatus  is  a  crossed  supersonic  molecular  beam 
device.  This  device  is  housed  under  a  high  vacuum  within  an  aluminum 
cylinder  four  feet  in  diameter  and  two  feet  high.  This  cylinder  is 
divided  into  three  distinct  chambers  as  shown  in  Figure  II-l.  The  first, 
and  the  largest,  is  the  main  chamber.  In  this  chamber  the  beams  intersect, 
the  fluorescence  is  detected,  and  both  beams  are  velocity  analyzed.  The 
second  and  third  chambers  are  virtually  identical  in  construction,  the 
only  differences  arising  in  the  additional  equipment  used  in  the  prepara- 
tion of  metastable  species.  Within  the  first  of  these  two  chambers,  the 
target  chamber,  is  a  nozzle  which  may  be  heated  by  a  resistive  element 
to  provide  variation  in  the  beam  velocity.  The  second,  and  last  of  the 
internal  chambers,  contains  a  nozzle  within  a  jacket  through  which  a 
coolant  liquid  may  be  passed  to  provide  for  velocity  variation  of  tnis 
beam.  Also  in  this  chamber  is  an  electron  gun  used  in  the  excitation 
of  the  ground  state  argon  atoms  to  their  metastable  states.  These  cham- 
bers are  each  differentially  pumped  by  oil  diffusion  pumps  with  typical 
pressures  attainable  listed  in  Table  II-l. 

The  gas-handling  system  has  been  extensively  described  by  Sanders 
[41]  and  Cutshall  [42]  and  so  will  not  be  detailed  in  this  paper.  The 
data  acquisition  system  and  the  calibration  system  for  the  target  beam 
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are  new  elements  to  the  experiment  and  as  such  will  be  extensively 
discussed. 

Table  II-l.  Typical  Vacuum  Pressures  in  Experimental  Apparatus 


Chamber 


Pressure  before  Experiment 


During  Experiment 


Main 

Metastable 

Target 


5  X  10"  torr 
1  X  10"^  torr 
1  X  10"  torr 


1  X  10"  torr 
5  X  10"^  torr 


5  X  10   torr 


B.  Gas-Handling  System 
The  purpose  of  the  gas-handling  system  is  to  provide  a  means  of 
admitting  high-purity  gases  into  the  vacuum  system  in  a  controlled 
manner.  By  means  of  a  series  of  valves,  shown  in  Figure  II-2,  the  gases 
are  reduced  in  pressure  from  several  hundred  atmospheres  in  a  standard 
high-pressure  cylinder  to  approximately  one  atmosphere  behind  the  nozzle 
orifice.  This  pressure  is  monitored  by  a  Wallace  and  Tiernan  model 
FA160  absolute  pressure  gauge  to  guarantee  beam  operation  under  constant 
conditions.  The  gas  then  expands  through  the  orifice  as  described  in 
Chapter  I  and  is  ultimately  pumped  out  of  the  system  and  vented  to  the 
atmosphere. 


C.  Metastable  Beam  Production 
The  production  of  the  metastable  beam  involves  conventional  use 
of  the  methods  of  supersonic  nozzle  beams  as  described  in  Chapter  I  but 
with  a  novel  source  for  exciting  the  metastable  states  in  the  beam  gas. 
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This  source  was  developed  by  Cutshall  [43]  for  previous  investigations 
of  the  Ar*/N2  reaction  system.  The  source  is  detailed  in  Figure  II-3, 
and  consists  of  a  tungsten  filament  cathode  located  to  the  side  of  the 
nozzle  cap  and  the  skimmer  opening.  The  electrons  are  boiled  off  the 
filament  by  a  15-ampere  current  passing  through  the  filament.  The 
electrons  are  accelerated  across  the  beam  at  right  angles  to  the  direc- 
tion of  the  beam  path  and  collected  by  a  nickel  plate  anode.  This  anode 
is  maintained  at  approximately  +30  volts  relative  to  the  filament  by  a 
constant  current  power  supply  and  regulates  the  discharge  current  at 
150  milliamperes.  This  source  is  capable  of  exciting  the  ground  state 
argon  atoms  into  a  variety  of  excited  states  as  well  as  ionizing  the 
atoms.  However,  because  of  the  length  of  the  flight  path  most  of  the 
excited  states  decay  prior  to  leaving  the  metastable  chamber,  and  the 
charged  particles  are  swept  out  of  the  beam  by  a  100-volt  potential 
applied  between  a  pair  of  sweep  plates  located  just  downstream  of  the 
skimmer  exit.  In  the  case  of  the  work  with  N^  as  the  target  gas,  the 
beam  then  proceeds  through  the  interaction  region,  is  detected  on  a 
Bendix  model  306  magnetic  electron  multiplier,  and  time-of-f light  analyzed 
to  determine  the  velocity  distribution  of  the  beam.  In  the  work  involving 
the  H2  target,  the  electron  multiplier  is  run  as  a  simple  surface  detec- 
tor from  which  electrons  are  ejected  by  the  metastable  atoms.  These 
electrons  are  then  detected  by  placing  a  10  -ohm  resistor  in  the  path 
to  ground  and  measuring  the  voltage  drop  across  the  resistor  by  an  elec- 
trometer. In  this  work  the  velocity  of  the  argon  metastable  atoms  was 
a  minor  contribution  to  the  relative  energy.  The  metastable  argon  beam 
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was  generated  at  the  same  source  conditions  used  in  the  room  temperature 
No  experiments,  and  the  velocity  distribution  was  assumed  to  be  the  same 
as  in  those  experiments. 

D.  Velocity  Analysis 

The  velocity  analysis  of  both  beams  is  carried  out  in  an  identical 
manner  with  only  the  particular  detection  method  varying.  For  the  pur- 
poses of  this  discussion,  the  specific  method  of  detection  is  irrelevant 
and  will  be  considered  as  a  black  box  from  which  emerges  a  signal  pro- 
portional to  the  number  of  particles  reaching  the  detector  per  second. 
It  should  be  mentioned  that  for  the  calculation  of  cross  sections,  the 
quadrupole  measures  a  number  density  while  the  electron  multiplier 
measures  an  intensity. 

The  beam  under  analysis  is  periodically  interrupted  by  a  slotted 
disk,  shown  in  Figure  1 1-4.  This  disk  is  driven  at  50  Hz  by  a  wide  band 
amplifier  which  is  in  turn  driven  by  an  oscillator.  This  disk  is  so 
constructed  as  to  allow  four  pulses  to  pass  the  chopper  in  every  cycle, 
two  long  pulses  and  two  short  pulses.  The  long  pulses  are  used  in  the 
photon  counting  portion  of  the  experiment  and  are  not  used  in  the  velocity 
analysis.  Only  the  short  pulses  are  used  to  prevent  a  smearing  out  of 
the  distribution  due  to  the  large  bandwidth  introduced  by  the  long 
pulse. 

As  the  beam  is  chopped,  so  too  is  a  light  path  between  a  light 
emitting  diode  (Texas  Instruments  TI-351)  and  a  phototransistor  (Texas 
Instruments  LS-400).  As  in  the  case  of  the  beam,  the  light  signal  is 
also  divided  into  two  separate  pulses,  an  optical  long  and  an  optical 
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short  pulse.  These  pulses  are  separated  by  electronic  circuitry  and 
used  independently.  The  relationship  between  the  optical  pulses  and 
the  beam  pulses  is  shown  in  Figure  II-5.  The  signals  are  then  fed  into 
a  Princeton  Applied  Research  model  121  phase-locked  amplifier  with  the 
optical  short  signal  being  used  as  the  reference  signal.  Signal  averaging 
is  accomplished  in  a  Princeton  Applied  Research  model  TDH-8  wave-form 
eductor  and  then  displayed  on  an  oscilloscope  triggered  by  a  delayed 
optical  short  pulse.  This  allows  the  intensity  to  be  measured  on  the 
lock-in  amplifier  during  the  velocity  analysis  portion  of  the  experiment. 

E.  Optical  System 
The  photon  collection  system  is  composed  of  a  mirror,  two  lenses, 
a  monochromator,  a  third  lens,  and  a  photomultiplier  positioned  as  shown 
in  Figure  II-6.  The  mirror  is  a  20-mm  (Rolyn  Optics  61.2200)  focal-length 
mirror  with  an  aluminum  overlay  coating.  This  mirror  is  located  one 
focal  length  below  the  interaction  region.  This  mirror  collects  all 
photons  emitted  in  the  downward  direction  and  reflects  them  back  through 
the  interaction  region.  Located  approximately  11  cm  above  the  inter- 
action region  is  a  100-mm  (Rolyn  Optics  11.2300)  focal-length  lens  which 
collimates  the  photons  into  parallel  rays  and  directs  them  up  the  2.5  in. 
o.d.  aluminum  tube  which  serves  as  a  mounting  structure  for  the  optical 
system.  Before  the  photons  reach  the  monochromator  the  photons  pass 
through  a  second  lens  with  a  focal  length  of  200  mm  (Rolyn  Optics  11.2650), 
which  serves  to  focus  the  light  on  the  entrance  slit  to  the  monochromator. 
The  photons  then  leave  the  vacuum  system  through  a  1.37  in.  diameter 
window  (Ceramaseal  #94764901-1).  The  monochramator  is  a  Spex  1670 
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Figure  II-6.  Photon  Collection  System 
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Minimate  mo no chroma tor  with  a  grating  ruled  at  1200  grooves  per  inch 
and  blazed  at  300  nm.  This  spectrometer  has  replaceable  slits.  A  1.0- 
nm  resolution  was  used  for  the  Np  work  and  10.0-nm  resolution  for  the 
Hp  experiments.  The  light  from  the  exit  slit  is  then  defocused  by 
lens  3,  a  25-mm  (Rolyn  Optics  11.2050)  lens.  This  defocusing  allows 
the  fullest  use  of  the  photocathode  surface.  All  lenses,  as  well  as  the 
window,  are  of  Herasil  grade  quartz.  The  detector  is  an  EMI  photomulti- 
plier  with  an  S-20  extended  photocathode.  The  photomultiplier  is  cooled 
to  -25°C  by  a  Products  for  Research  104-TE-RF  thermoelectric  refrigera- 
tor. This  reduces  the  dark  current  which  would  otherwise  be  the  major 
source  of  noise  in  the  system. 

F.  Data  Collection  System 

After  a  photon  has  been  detected  by  the  photomultiplier,  the 
pulse  must  be  analyzed  in  terms  of  the  chopped  beam  phases.  This  is 
accomplished  by  the  digital  electronics  referred  to  in  Figure  II-8  as 
the  MCA/Stepping  Motor  Interface.  This  figure  also  traces  the  path  of 
the  photon  signal  from  the  photomultiplier  to  the  storage  or  rejection 
of  the  pulse.  The  pulse  is  first  amplified  by  a  Tennelec  TC-145  pre- 
amplifier, further  amplified  by  a  Tennelac  TC-202  linear  amplifier, 
filtered  and  converted  to  a  square  pulse  by  a  Tennelac  TC-440  single 
channel  analyzer,  and  finally  put  out  to  the  interface  circuitry  where 
the  signal  is  compared  to  the  signal  generated  by  the  chopper  mechanism 
as  the  target  beam  is  periodically  interrupted. 

This  signal,  the  optical  long,  is  first  delayed  in  a  simple  delay 
circuit,  shown  in  Figure  II-7: 
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Figure  II-7.  Pulse  Delay  Circuit 


This  circuit  takes  advantage  of  the  intrinsic  delay  times  through  the 
exclusive-or  gates  to  provide  an  output  pulse  of  short  duration  (50  nsec) 
on  each  change  of  state  of  the  input  signal  corresponding  to  the  opening 
and  closing  of  the  chopper.  This  signal  then  triggers  a  monostable 
vibrator  (TTL-74123)  which  outputs  a  pulse  after  a  time  delay  controlled 
by  a  ten-turn  potentiometer  adjustable  from  the  front  panel.  This  pulse 
then  triggers  a  D-type  flip-flop  (TTL-7474)  which  outputs  a  signal 
dependent  upon  the  current  state  of  the  chopper. 

This  circuit  puts  out  two  complementary  signals  which  are 
delayed  by  a  given  time  interval  from  the  opening  of  the  chopper  to 
allow  the  target  beam  to  reach  a  steady  state  within  the  interaction 
region  prior  to  beginning  data  counting.  The  signal  which  takes  on  the 
same  state  as  the  original  optical  long  then  serves  a  triple  function. 
This  signal  feeds  the  interface  to  provide  a  signal  which  is  used  to 
synchronize  all  of  the  timing  within  the  interface,  triggers  an 
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Elscint  CBC-N-1  crystal  clock  which  opens  a  window  in  the  interface 
during  a  set  time  interval  in  which  the  beam  is  actually  in  the 
interaction  region,  and  finally  triggers  a  Tennelec  TC-551  preset 
scalar  which  monitors  the  stepping  function  of  the  system. 

The  complementary  signal  is  fed  directly  to  the  data  storage 
device,  a  Nuclear  Data  series  1100  multichannel  analyzer  (MCA),  as  a 
reference  signal  determining  whether  a  given  photon  is  to  be  added  into 
the  memory  as  a  signal  +  noise  pulse  or  subtracted  from  memory  as  a 
noise  pulse. 

The  remaining  input  into  this  interface  is  the  output  of  a  free- 
running  BNC-8010  pulse  generator  which  controls  the  stepping  rate  of 
the  Slo-Syn  translator  and  stepping  motor. 

The  outputs  from  this  interface  are  the  photon  pulses  which 
arrive  in  the  predefined  time  intervals  (DATA  OUT),  a  pulse  which  tells 
the  MCA  to  advance  a  channel  (MCS  XTL  CLK) ,  an  output  which  keeps 
track  of  the  number  of  channels  stepped  to  allow  for  reset  of  the  sys- 
tem (SCALER),  a  pulse  train  telling  the  stepping  system  to  act  (UP/DN), 
and  a  pulse  which  informs  a  second  MCA  (Tracor  Northern  TN-1710)  which 
was  used  in  some  experiments  to  gather  signal-to-noise  data.  The 
complete  schematic  diagram  for  the  interface,  as  well  as  the  timing 
diagrams  for  the  input  and  output  signals,  are  included  in  Appendix  I. 
After  the  experiment  has  been  concluded  the  output  from  the  Nuclear 
Data  MCA  may  be  read  out  on  either  an  oscilloscope  or  as  hard  copy  on 
a  Teletype  ASR-33.  These  signal  trains  are  shown  in  Figure  II. 8. 
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G.  Mass  Spectrometer  Calibration  System 
As  seen  in  Chapter  I,  a  capillary  array  molecular  beam  is 
primarily  dependent  only  upon  the  pressure  of  the  gas  behind  the 
array.  This  guarantees,  if  the  applied  pressure  is  accurately  known, 
that  the  beam  will  be  reproducible  on  a  day-to-day  basis.  This  prop- 
erty makes  the  capillary  array  source  ideal  for  use  as  a  calibration 
source  for  the  quadrupole  mass  spectrometer  used  in  these  experiments. 

The  source,  shown  in  Figure  II-9,  is  composed  of  three  major 
components,  the  gas  introduction  lines,  the  pressure  measuring  device, 
and  the  capillary  array  itself.  The  gas  introduction  system  is  simply 
a  1/4  in.  copper  line  tapped  into  the  target  gas  line  and  terminated 
by  a  leak  valve.  This  leak  valve  is  used  to  control  the  pressure  be- 
hind the  capillary  array.  For  the  H  /Ar*  cross  section  measurements 
this  pressure  was  regulated  at  0.500  torr  relative  to  the  main  chamber 
pressure  ('v^  1  x  10"^  torr).  This  pressure  was  measured  by  an  MKS 
Instruments,  Inc.  Baratron  Model  90  capacitive  manometer.  One  end  of 
the  device  was  attached  to  the  high  pressure  region  behind  the  cap- 
illary array  and  the  other  end  opened  into  the  main  vacuum  chamber. 
These  two  sides  were  connected  by  a  valve  which  allowed  for  rapid 
pump  out  of  the  high  pressure  region. 

The  beam  source  itself  was  located  in  a  small  block  of  aluminum 
and  contained  the  capillary  array,  a  chopping  mechanism,  and  an  LED- 
phototransisotr  system  identical  to  that  used  with  the  supersonic 
sources.  The  capillary  array  maintained  the  pressure  difference  be- 
tween the  source  and  the  main  chamber.  It  also  provided  the 
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conimation  required  for  high  beam  intensities.  As  with  the  super- 
sonic nozzles,  this  beam  was  chopped  providing  short  pulses  of  gas 
which  were  detected  by  the  quadrupole  mass  spectrometer.  The  design 
of  the  chopper  mechanism  and  associated  electronics  was  such  that  the 
detection  procedures  were  identical  to  those  discussed  in  the  time-of- 
f light  measurements. 


CHAPTER  III 
EXPERIMENTAL  PROCEDURE 


A.  Optical  Considerations  for  Ar*/N2 
The  flowing  afterglow  experiments  of  Setser  and  Stedman  [44]  and 
Setser,  Stedman  and  Coxon  [45]  have  shown  three  emission  bands  for  the 
Ar*/Np  system.  These  bands  correspond  to  the  transitions 

c  -^n  -i-  B  ^n^ 

"      9       2nd  Positive 

B  "^n  -^A  ^Z'*'       1st  Positive 
9      u 

A  ^Z"^  -^X  ^E"^       Vegard-Kaplan 
u      g 

Figure  III-l  shows  the  relevant  potential  energy  curves  for  N  and  N  . 
Interference  with  the  desired  C^^  transitions  are  not  important  in  this 
work  as  the  A^X  transition  is  metastable,  transitions  occurring  several 
milliseconds  downstream.  This  time  lag  causes  all  A->X  fluorescence  to 
occur  outside  of  the  detection  region.  The  B^A  transition  is  short 
lived;  however,  this  radiation  occurs  in  the  infrared  region  of  the  spec- 
trum and  does  not  pass  the  monochromator. 

Of  primary  importance  to  this  experiment  are  the  0-0  band  and 
the  1-0  band  of  the  C^B  electronic  transition.  These  are  found  to  lie 
in  the  near  ultraviolet  in  the  ranges  330-339  nm  and  310-318  nm. 
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Figure  III-l.  Np  Potential  Energy  Curves 
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respectively.  These  regions  were  repetitively  scanned,  using  1  nm 
spectrometer  resolution,  until  the  peak  channel  contained  in  excess  of 
2,000  counts.  This  value  gave  a  modest  signal-to-noise  ratio,  approxi- 
mately 2:1,  within  a  reasonable  data  collection  time. 

B.  Optical  Considerations  for  Ar*/H2 

James  and  Coolidge  [46]  have  calculated  the  anticipated  H2  con- 
tinuum spectra.  For  the  v'  =  0  vibrational  state,  the  emission  is 
expected  to  lie  between  220  nm  and  approximately  500  nm,  with  a  maxi- 
mum emission  at  260  nm.  The  appropriate  potential  energy  curves  for 
the  H2  molecule  [47]  are  given  in  Figure  III-2. 

An  alternative  source  of  photons  is  the  continuum  emission  from 
an  excited  bound  state  of  an  ArH  molecule  to  the  unbound  ground  state. 
Figure  I I 1-3  shows  the  potential  energy  curves  for  this  molecule  as 
calculated  by  Olson  [48].  It  can  be  seen  that  the  emission  would  cover 
approximately  the  same  region  of  the  spectrum,  although  a  lower  energy 
maximum  may  be  expected. 

Regardless  of  the  product,  the  region  to  be  explored  ranged 
from  220  nm  to  500  nm.  This  region  was  scanned  only  one  time,  but  with 
a  10  nm  bandwidth  and  a  10  second  dwell  time  per  step. 

C.  Experimental  Procedure 
The  first  step  in  the  experimental  procedure  was  the  evacuation 
of  the  vacuum  chambers.  The  entire  system  was  rough  pumped  by  the 
main  chamber  backing  pump  through  the  master  valve  and  three  Veeco 
valves  located  directly  beneath  each  internal  chamber.  During  this 
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time  the  diffusion  pumps  were  isolated  from  their  respective  backing 
pumps.  After  the  pressure  in  the  chambers  fell  below  750  microns,  the 
refrigerator  to  the  main  chamber  baffle  was  turned  on.  After  at  least 
15  minutes  and  after  the  pressure  in  the  system  had  fallen  below  100 
microns,  the  cooling  water  to  the  three  diffusion  pumps  and  the  baffle 
to  the  metastable  chamber  was  turned  on  and  the  backing  lines  to  the 
metastable  chamber  and  target  chamber  were  evacuated.  At  this  point 
the  diffusion  pumps  were  connected  with  the  backing  lines  by  opening  the 
isolation  valves.  The  diffusion  pumps  were  turned  on  and  the  Veeco 
valves  and  master  valve  closed.  The  system  was  allowed  to  pump, 
usually  overnight,  until  the  static  pressures  noted  in  Chapter  II  were 
attained. 

After  the  system  had  been  evacuated,  it  was  necessary  to 
establish  the  beams,  prepare  the  electronics,  and  ready  the  photon 
counting  system.  Adequate  warmup  time  (several  hours)  for  the  electronics 
to  stabilize  was  necessary.  The  critical  electronics  include  the  power 
supplies  for  the  metastable  discharge,  the  oscillators  driving  the 
choppers,  the  high-voltage  power  supplies  for  the  quadrupole  mass 
spectrometer,  electron  multiplier,  and  the  photomultipler,  the  electrome- 
ter, lock-in  amplifier,  and  wave-form  eductor.  All  other  electronic 
elements  did  not  require  extended  warmups,  although  they  were  normally 
turned  on  at  the  same  time  as  the  critical  devices. 

The  establishment  of  the  beams  required  the  same  procedure  for 
either  beam,  although  generation  of  the  metastable  beam  required  the 
additional  step  of  generating  the  electrical  discharge  between  the  nozzle 
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and  skimmer.  The  first  step  in  setting  up  the  beams  was  the  evacuation 
of  the  nozzle  lines.  After  the  line  pressure  fell  below  50  microns  the 
vacuum  pumps  were  isolated  from  the  system.  With  the  valve  connecting 
the  nozzle  to  the  gas  line  open  and  the  leak  valve  and  leak-valve  bypass 
valve  closed  the  gas  was  allowed  to  enter  the  line  at  20  psig  through  a 
regulator  attached  to  a  high-pressure  gas  cylinder.  The  leak  valve  was 
then  slowly  opened,  allowing  the  gas  to  flow  to  the  nozzle.  The  pressure 
behind  the  nozzle  was  stabilized  at  15  psia  by  adjusting  the  leak  valve 
setting.  If  the  experiment  required  the  heating  or  cooling  of  the 
nozzle,  the  procedure  was  begun  at  this  time.  The  target  beam  was  now 
allowed  to  stabilize,  a  process  requiring  about  two  hours,  prior  to  the 

beginning  of  data  collection.  As  the  pressure  in  the  metastable  chamber 

-4 
rose  to  approximately  1  x  10   torr,  the  electron  gun  was  turned  on  and 

the  power  supplies  adjusted  to  provide  150  milliamperes  discharge  cur- 
rent. The  metastable  beam  was  then  allowed  several  hours  to  stabilize. 
Stabilization  of  the  beams  was  detected  by  a  steady  backing  pressure  to 
the  nozzle,  a  steady  output  from  the  associated  detection  device,  and 
most  importantly  by  constant  measurements  of  the  beam  velocity. 

While  the  beams  were  stabilizing,  the  photon  collection  system 
was  prepared.  This  required  the  refrigeration  of  the  photomultiplier, 
the  heating  of  the  photomultiplier  window  to  prevent  condensation  on 
the  window,  the  setting  of  the  monochromator  to  its  initial  position, 
and  the  clearing  of  the  memory  of  the  MCA.  Just  prior  to  the  photon 
counting  portion  of  the  experiment  liquid  nitrogen  was  added  to  a  cryo- 
pump,  thereby  reducing  the  background  due  to  condensibles  in  the  vacuum 
system. 
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Photon  counting  was  initiated  by  opening  the  shutter  to  the 
photomultiplier,  chopping  the  target  beam,  and  pressing  the  start  but- 
ton on  the  MCA/Stepping  Motor  Interface  panel.  After  sufficient  data 
were  collected,  the  process  was  terminated  by  pressing  the  stop  button 
on  the  same  panel.  The  photon-counting  portion  of  the  experiment  then 
ceased  to  be  operative  at  the  end  of  the  current  scan. 

At  the  end  of  a  run,  a  measurement  of  the  velocity  of  the  beams 
by  means  of  the  time-of-f light  technique  was  made.  The  short  burst  of 
particles  passing  through  the  small  slot  in  the  chopper  was  detached 
by  the  quadrupole  in  the  case  of  the  target  beam,  or  the  electron  mul- 
tiplier for  the  metastable  beam.  This  signal  was  amplified  by  the  lock- 
in  amplifier  and  signal  averaged  by  the  wave  form  eductor.  The  signal 
was  then  compared  for  timing  purposes  with  the  optical  short  signal  and 
a  simple  distance-over-time  calculation  provided  the  most  probable 
velocity  of  the  beam.  For  cross-section  measurements,  an  analog  signal 
from  the  lock-in  amplifier,  which  was  proportional  to  the  target  beam 
intensity  was  measured.  The  proportionality  constant,  which  depended 
on  daily  operating  conditions,  was  determined  by  producing  a  known 
pressure,  as  measured  by  the  capacitance  manometer,  behind  the  capillary 
array.  The  resulting  beam  was  then  monitored  in  the  same  fashion  as  the 
target  beam  had  been  monitored. 

Other  correction  factors  are  given  in  Figures  III-4  and  III-5. 
Figure  III-4  is  a  wavelength  calibration  of  the  spectrometer  system 
determined  using  a  mercury  vapor  lamp.  Figure  III-5  is  an  intensity 
calibration  based  on  a  comparison  with  a  deuterium  standard  lamp. 
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CHAPTER  IV 
DATA  ANALYSIS  AND  RESULTS 


A.  N2  Band  Profiles  and  Rotational  Distributions 

The  profiles  of  the  N2(C^n^  (v'  =  0))  to  N2(B^n  (v"  =  0)  band 
were  measured  at  the  relative  collision  energies  of  0.076  eV,  0.080  eV, 
and  0.161  eV.  These  profiles  are  given  in  Figure  IV-1  as  the  vertical 
slash  marks,  the  size  of  these  marks  being  determined  by  the  estimated 

errors  on  the  particular  data  point.  Figure  IV-2  shows  the  profiles  for 

3  3 

the  N2(C  HjjCv'  =  1)  to  N2(B  n  (v"  =  0))  transitions  obtained  at  the  same 

relative  energies  as  the  0-0  bands. 

The  smooth  curves  through  the  data  points  on  the  two  sets  of 
spectra  represent  the  best  fit  to  the  data  based  on  known  spectroscopic 
constants  and  various  trial  rotational  distributions.  The  procedure  for 
obtaining  these  fits  was  to  assume  a  particular  rotational  distribution, 
feed  this  into  a  computer  program  [49]  along  with  the  spectroscopic 
constants  for  the  relevant  states,  broaden  the  spectra  to  allow  for  the 
spectrometer  bandwidth,  and  finally  generate  a  least-squares  fit  to  the 
original  data.  This  procedure  was  then  repeated  with  a  new  trial  distri- 
bution of  rotational  states  until  a  minimum  was  obtained  in  the  least- 
squares  fit.  The  resultant  spectrum  was  then  compared  visually  to  the 
data  and  a  decision  was  made  as  to  the  advisability  of  attempting  the 
fitting  procedure  with  a  different  type  of  distribution.  A  listing  of 
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the  computer  program  used  to  generate  these  spectra  is  included  in 
Appendix  II. 

For  the  0-0  band  it  was  found  that  the  data  could  be  adequately 
fit  by  a  thermal,  i.e.,  Boltzmann,  distribution.  As  this  entailed  only 
the  fitting  of  a  "temperature"  to  the  data  the  routines  for  the  fitting 
were  quite  easily  accomplished.  The  results  of  this  procedure  are  given 
in  Figure  IV-3. 

For  the  1-0  band  a  simple  thermal  distribution  was  not  found  to 
adequately  fit  the  data  and  more  complex  procedures  were  required.  The 
most  successful  procedure  involved  calculating  the  spectra  using  the 
first  several  Chebyshev  polynomials  [50]  with  arbitrary  coefficients  to 
represent  the  rotational  distribution,  and  then  performing  a  correlation 
of  the  results  to  determine  the  best  coefficients  to  be  used.  Upon 
visual  inspection  a  second  distribution  suggested  itself.  This  is  the 
so-called  "Golden  Rule"  distribution  discussed  in  the  Introduction, 
Equations  (30)  and  (31).  These  distributions  are  both  plotted  in  Figure 
IV-4. 

Summarizing  these  results  given  in  Figure  IV-3,  the  0-0  band  is 
in  all  cases  characterized  by  a  rotational  "temperature."  These  "tem- 
peratures" are  2400  K,  1600  K,  and  1700  K,  and  have  an  estimated  error  of 
100  K.  The  1-0  band  may  be  represented  by  either  a  series  of  Hermite 
polynomials  or  by  a  "Golden  Rule"  distribution.  In  the  case  of  the 
Chebyshev  polynomials  the  distribution  is  given  by 

P(J)  =AxTq+BxT^+CxT2 
where  the  Chebyshev  polynomials  on  the  interval  0-1  are 
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The  series  was  cut  off  based  on  the  statistics  which  indicated  that 
no  further  information  was  gained  with  additional  terms. 

B.  Hj,  +  Ar*  Cross  Sections 

and  Spectral  Distributions 
The  behavior  of  the  cross  section  as  a  function  of  the  relative 
translational  energy  of  the  collision  partners  is  quite  typical  of  any 
endoergic  process  in  that  there  is  no  contribution  to  the  cross  section 
until  a  certain  minimum  energy  is  obtained,  after  which  there  occurs  a 
rapid  rise  to  a  maximum  and  then  a  slow  tapering  off.  Figure  IV-5  is 
the  experimentally  determined  energy  dependence  of  the  cross  section 
for  the  process  involving  the  collision  of  H2(X'ag)  +  Ar*(  P2  q).  The 
error  bars  for  the  energy  determination  represent  the  full  width  half- 
maximum  distribution  of  the  beam  velocities  and  the  error  bars  on  the 
value  of  the  cross  section  represent  95%  confidence  intervals.  All  data 
points  have  been  calculated  according  to  Equation  (19)  in  Chapter  I. 

The  major  characteristic  of  this  curve  is  the  rapid  rise  in  the 
cross  section  beginning  at  the  onset  at  0.08  eV  and  the  apparent  maximum 
at  0.170  eV.  Because  of  limitations  in  the  design  of  the  experiment, 
higher  energies  were  not  attainable  and  the  behavior  at  higher  energies 
could  not  be  ascertained. 
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Figure  IV-5.  Energy  Dependent  Cross  Section  for  Ar*  +  H2  Collisions 
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Figures  IV-6  and  IV-7  are  spectral  distributions  of  the  process 
taken  at  the  lowest  (0.067  eV)  and  highest  (0.165  eV)  relative  energies, 
respectively.  These  distributions  have  been  smoothed  by  a  three-point 
smoothing  routine  [51]  and  corrected  for  spectral  response  of  the 
photon-gathering  system.  An  estimate  of  the  errors  on  the  individual 
data  points  near  the  peak  of  the  distributions  was  obtained  by  position- 
ing the  spectrometer  on  the  peak  and  running  the  experiment  without 
stepping  the  spectrometer.  These  data  were  then  run  through  the  same 
smoothing  procedure  and  the  errors  estimated  from  the  variation  observed 
in  these  data.  It  was  found  that  the  errors  on  the  data  points  are 
approximately  the  same  as  the  apparent  oscillations  in  the  distribution. 

The  results  of  the  photon  distribution  experiment  indicate  that 
the  maximum,  regardless  of  collision  energy,  lies  at  approximately  290  t 
5  nm  and  that  the  distributions  of  photons  are  nearly  identical  over  all 
observed  wavelengths  independent  of  collision  energy. 
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CHAPTER  V 
DISCUSSION 

A.  Collisions  Involving  Ar*/N2 

3 
The  rotational  distributions  observed  for  the  N^CC  IT,  v'  =0) 

product  of  the  Ar*/N2  collision  system  may  be  represented  well  by 
Boltzmann  distributions  over  the  range  of  collision  energies  of  0.078  eV 
to  0.161  eV.  The  Boltzmann  temperatures  for  these  distributions  vary 
from  1600  K  to  2200  K.  For  these  collisions  there  is  sufficient  energy 
to  populate  rotational  levels  up  to  the  45th  to  50th  quantum  level.  The 
distributions,  as  shown  in  Figure  IV-3,  rise  rapidly  to  maximum  popula- 
tions in  the  range  J  =  18-20,  and  then  slowly  taper  off  until  the 
maximum  allowed  J-value  is  reached. 

In  contrast  to  the  v'  =  0  levels,  the  N2(C  11,  v'  =  1)  product 
shows  a  marked  difference  in  the  distribution.  In  fact  the  spectra 
could  not  be  fit  by  a  Boltzmann-type  distribution.  Instead,  the  spectra 
were  found  to  be  reproduced  within  the  experimental  error  by  the 
"Golden  Rule"  distributions  shown  in  Figure  IV-4.  These  distributions, 
because  of  the  additional  energy  required  to  enter  the  v'  =  1  vibrational 
level,  do  not  have  sufficient  energy  to  populate  higher  rotational  levels 
and  the  distributions  cut  off  at  approximately  J  =  35.  The  character- 
istics of  these  distributions  are  a  more  gentle  rise  to  the  maximum 
population  levels  of  J  =  25-30  and  then  slowly  begin  to  taper  off.  Be- 
cause the  cut-off  occurs  so  quickly  after  the  maximum  is  reached,  the 
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relative  population  levels  are  still  quite  high  at  the  maximum  allowable 
value  for  J. 

No  good  explanation  as  to  why  the  v'  =0  vibrational  level 
appears  to  be  represented  by  a  Boltzmann  distribution  is  currently  avail- 
able. Although  such  a  distribution  would  be  expected  in  the  case  in 
which  the  product  was  in  thermal  equilibrium,  in  the  single-collision 
processes  examined  in  the  molecular  beam  experiment  no  such  equilibra- 
tion occurs. 

B.  Collisions  Involving  Ar*/H2 

The  photon  energy  distribution  data  and  the  colli  si onal  energy 
dependence  of  the  cross  section  for  the  Ar*/\\^   system  do  not  unambigu- 
ously define  the  process  being  observed.  Two  possible  processes  may  be 
responsible,  either  singly  or  in  tandem,  for  the  observed  results. 
These  processes  are  the  transfer  of  electronic  energy  to  the  Hp  diatomic 

q 

with  the  observed  continuum  radiation  corresponding  to  the  H2(a  n  , 

q 

v'  =  0)  — »-  H2(b  n  )  transition  or  an  exchange  process  resulting  in  a 
bound  excited  state  of  ArH*  which  then  fluoresces  to  the  unbound  ground 
state,  again  with  a  broad-band  emission  of  photons.  Each  of  these 
processes  will  be  discussed  with  respect  to  the  consequences  of  the 
particular  assumption. 

B.l  Radiation  from  H2(a  ^n^,  v'  =  0) 

The  process  of  exciting  the  ground  state  molecule  to  the  excited 

q 

"a"  state  requires  11.788  eV  of  energy.  The  Ar*(  Pq)  metastable  stores 
11.723  eV  resulting  in  a  deficit  of  0.065  eV  which  must  be  the  minimum 
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amount  of  energy  provided  by  the  translational  motion  of  the  system  for 
the  process  to  occur.  The  Ar*(  P-)  metastable  provides  only  11.548  eV 
to  the  system  with  a  resulting  deficit  of  0.240  eV.  As  the  maximum 

energy  available  from  translational  motion  in  the  experiment  is  0.161  eV, 

3 
the  Pp  component  cannot  contribute  to  the  process.  Also  the  v'  =  1 

vibrational  level  of  the  "a"  state  lies  approximately  0.243  eV  above  the 

v'  =  0  level,  making  this  level  energetically  unattainable  by  any  species 

3 
in  the  system.  Although  the  H„(c  n  ,  v'  =  0)  is  attainable  energetically, 

^    y 

the  state  is  metastable  and  does  not  decay  in  the  detection  zone. 

From  the  above  results,  it  is  possible  to  predict  the  behavior 
of  the  energy-dependent  cross  sections  for  the  ArVHp  system.  The  cross 
section  should  be  zero  until  the  relative  translational  energy  reaches 
0.065  eV,  after  which  the  cross  section  should  rise  rapidly  to  a 
maximum  and  then  slowly  taper  off.  This  does  assume  that  no  activation 
energy  is  required,  in  which  case  the  onset  would  occur  at  a  still  higher 
relative  energy.  Examination  of  Figure  IV-5  reveals  this  type  of  a 
behavior  with  the  onset  occurring  at  approximately  0.080  eV  implying 
an  activation  energy  of  0.015  eV. 

Serious  questions  about  this  interpretation  arise  when  the 
spectral  distributions  are  examined.  As  seen  in  Figure  IV-6  the  ob- 
served distributions  are  significantly  different  than  those  predicted 
by  James  and  Coolidge  [51].  The  potential  energy  curves  calculated  in 
this  work  compare  quite  well  with  the  later  work  of  Kolos  and  Wolniewicz 
[52],  therefore  the  theoretical  prediction  must  be  considered  reliable. 
The  observed  spectral  distributions  might  arise  from  potential  energy 
curves  perturbed  by  the  presence  of  the  heavy  argon  atom.  This 
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hypothesis  runs  into  trouble  when  the  observed  distributions  at  differ- 
ent relative  energies  are  compared.  At  the  higher  relative  energies, 
it  would  be  expected  that  the  argon  atom  would  be  further  away  from 
the  hydrogen  molecule  at  the  time  of  the  radiation.  The  perturbation 
would  thus  be  smaller  and  the  distribution  of  photons  would  be  differ- 
ent. This  does  not  agree  with  the  spectra  observed  in  Figure  IV-6  and 
Figure  IV-7,  which  show  nearly  identical  distributions. 

One  further  piece  of  evidence  damaging  to  this  model  is  found 
in  the  work  of  Feldstein  [53].  He  found  that  the  onset  for  radiation 
in  the  Ar*/D2  system  lies  at  0.075  eV,  0.023  eV  below  the  expected  onset 
of  0.098  eV,  a  situation  in  which  the  "a"  state  of  the  D^  molecule 
cannot  possibly  be  populated. 

B.2  Radiation  from  ArH* 

The  ArH*  molecule  has  been  observed  in  the  visible  region  of 
the  spectrum  by  Johns  [54].  The  potential  energy  curves  have  been 
calculated  by  Olson  [55]  using  SCF-CI  methods  (see  Figure  IV-8).  Using 
the  potential  energy  parameters  derived  from  these  sources,  it  is 

possible  to  calculate  the  endoergicity  for  the  formation  of  the 

2  3 

ArH*(A  E,  v'  =  0)  molecule.  For  the  Pq  metastable  the  endoergicity 

3 
is  -0.72  eV,  and  for  the  Pp  metastable  the  endoergicity  is  -0.55  eV. 

In  all  cases  the  process  should  occur  at  all  relative  energies,  if  no 

activation  energies  are  involved.  The  cross  section  data  previously 

mentioned  can  then  be  explained  by  an  activation  energy  on  the  order  of 

0.70  eV.  By  being  exoergic,  the  results  of  Feldstein  for  the  ArVDp 

system  are  no  longer  in  conflict  with  the  anticipated  results.  A 
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comparison  of  the  energy  differences  between  the  excited  bound  states 
and  the  lower  energy  unbound  states  at  the  equilibrium  positions  of 
the  bound  states 

H2(a-^b)  rg  =  0.9888  A       AE  =  4.02  eV 

ArH(A-^X)  rg  =  1.2686  A       AE  =  3.58  eV 

indicate  that  the  spectral  distribution  would  be  expected  to  peak  at  a 
higher  wavelength  for  the  ArH  transitions.  The  exact  position  of  this 
maximum  cannot  be  predicted  without  the  calculation  of  the  overlap 
integrals  between  the  ArH(A)  and  the  unbound  ArH(X)  states. 

B.3  Conclusion 

The  current  data  tend  to  support  the  premise  that  the  process 
being  observed  is 

Ar*(^P2  q)  +  H2(X  ^Zp  -.  ArH*(A  h)   +  Ar(^S) 

although  the  production  of  H2{a  E^^)  cannot  be  ruled  out  as  contributing 
to  the  total  radiation  observed.  This  conclusion  may  be  further  tested 
experimentally  by  observing  the  process  in  which  the  H^  molecule  is 
excited  to  the  "a"  state  by  electron  impact  and  comparing  the  observed 
fluorescence  from  this  process  with  that  observed  in  this  study.  A 
calculation  of  the  three-body  potential  energy  surfaces  for  the  (Ar/Hp)* 
system  may  also  yield  information  making  the  interpretation  a  more 
definite  one. 

The  surfaces  of  interest  would  include  the  Ar*-H-H  surface.  This 
surface  represents  the  incoming  channels  of  the  collision  process  on 
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which  the  two  hydrogen  atoms  would  be  closely  associated  with  one 

another  and  the  Ar*  atom  would  be  approaching  the  pair.  This  surface 

then  must  cross  the  Ar-H(2s)*-H  surface  as  this  surface  contains  the         ''* 

two  primary  exit  channels  for  the  process.  Along  one  path  the  three-body 

collision  complex  would  exit  as  a  normal  Ar  atom  and  an  excited  hydrogen 

pair,  which  would  have  a  normal  hydrogen  atom  and  an  excited  (2s)  hydro-      ^ 

gen  atom  in  its  dissociation  limit.  The  other  path  would  have  as  its 

dissociate  limit  a  ground  state  ArH. 
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INTEGRATED  CIRCUIT  CHIP  LOCATION  AND  IDENTIFICATION  CHART 


CHIP 
#* 

INTEGRATED  CIRCUIT 
#  (TTL) 

FUNCTION 

1 

7405 

Hex  Inverter  (Open  Collector) 

2 

7404 

Hex  Inverter 

3 

7410 

3-Input  NAND  Gate 

4 

7400 

Quad  2-Input  NAND  Gate 

5 

7450 

And-Or- Invert  Gate 

6 

7420 

Dual  4- Input  NAND  Gate 

7 

7402 

Quad  2-Input  NDR  Gate 

8 

7474 

Dual  D  Flip-Flop 

9 

7404 

Hex  Inverter 

10 

7474 

Dual  D  Flip-Flop 

n 

7405 

Hex  Inverter  (Open  Collector) 

12 

7474 

Dual  D  Flip-Flop 

13 

7410 

3-Input  NAND  Gate 

14 

74122 

Monostable  Multivibrator 

15 

74123 

Monostable  Multivibrator 

16 

74190 

Decade  Up/Down  Counter 

17 

74190 

Decade  Up/Down  Counter 

18 

74190 

Decade  Up/Down  Counter 

19 

74190 

Decade  Up/Down  Counter 

*Chips  are  numbered  top  to  bottom,  left  to  right  when  the  board  is 
observed  with  the  edge  connector  up  and  the  chips  facing  toward  the 
observer. 
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MAIN 

Q********************************************************************* 

C       THIS  PROGRAM  WILL  CALCULATE  THE  SPECTRUM  FOR  A  TRANSITION  FROM  A 

C  3-PI-U  STATE  TO  A  3-PI-G  STATE  SUCH  AS  IN  THE  CASE  OF  THE 

C  N2(C-B)   TRANSITIONS.     FOR  TRANSITIONS   INVOLVING  OTHER  TYPES 

C  OF  SYMMETRY  STATES,  THE  HONL-LONDON  FACTORS  MUST  BE  ALTERED. 

C  THIS  SECTION  OF  THE  PROGRAM  IS  MARKED  IN  THE  BODY  OF  THE  TEXT. 

C  THE  PROGRAM  FIRST  CALCULATES  A  STICK  SPECTRUM  FROM  THE  SPEC- 

C  TROSCOPIC  DATA  SUPPLIED  AND  GIVEN  ROTATIONAL  DISTRIBUTIONS. 

C  THE  RESULTS  ARE  THEN  CORRECTED  FOR  SPECTROMETER  SENSITIVITY 

C  AND  SPECTRAL  RESOLUTION.     THE  RESULTS  ARE  THEN  PRINTED  OUT 

C  AND  IF  DESIRED  PLOTTED  USING  THE  GOULD  PLOT  ROUTINES. 

C 

C       REFERENCE:     THE  SPECTRA  OF  DIATOMIC  MOLECULES  BY  G.   HERZBERG. 

C 

C       THE  DATA  IS  INPUT  TO  THE  PROGRAM  IN  THE  FORM  OF  A  SERIES  OF  NAME- 

C  LIST  READS.     THE  NAMELIST  NAME,   DATA  NAME,  AND  A  DESCRIPTION 

C  OF  THE  DATUM  INCLUDING  THE  DEFAULT  VALUES  ARE  GIVEN  IN  THE 

C  TABLE  BELOW. 

C 

C       NAMELIST  VARIABLE               DESCRIPTION 

C 

C       UPPER OMEGA THIS   IS  AN  ARRAY  OF  DIMENSION  12, 

C  AND  CONTAINS  THE  FIRST  12  VIBRA- 

C  TIONAL  CONSTANTS  AS  DETERMINED  BY 

C  SPECTROSCOPIC  ANALYSIS.     THE  CON- 

C  STANTS  ARE  TO  THE  SERIES  OF  THE  FORM 

C  0MEGA(L)*(V+0.5)**L 

C  WHERE  L  IS  THE  SUMMATION  VARIABLE 

C  WHICH  RUNS  FROM  1   TO  12.     V   IS  THE 

C  VIBRATIONAL  QUANTUM  NUMBER  OF  THE 

C  STATE  UNDER  CONSIDERATION.    (DEFAULT=0) 

C 

C  ALPHA THIS  IS  AN  ARRAY  OF  DIMENSION  12, 

C  AND  CONTAINS  THE  FIRST  12  ROTATIONAL 

C  CONSTANTS  TO  THE  SERIES 

C  ALPHA(L)*(V+0.5)**L 

C  WHERE  L   IS  THE  SUMMATION  VARIABLE 

C  WHICH  RUNS  FROM  1   TO  12.     V   IS  THE 

C  VIBRATIONAL  CONSIDERATION.    (DEFAULT=0) 

C 

C  BE THIS  IS  THE  ROTATIONAL  CONSTANT  AT 

C  EQUILLIBRIUM  AND  IS  ASSOCIATED  WITH 

C  THE  ROTATIONAL  ENERGY  OF  THE  MOLECULE 

C  BE*J*(J+1) 

C  WHERE  J   IS  THE  ROTATIONAL  QUANTUM 

C  NUMBER.      (DEFAULT=0) 

C 

C  DE THIS  IS  THE  ROTATIONAL  CONSTANT  AT 

C  EQUILLIBRIUM  AND  IS  ASSOCIATED  WITH 


■4', 
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C  THE  SECOND  ORDER  CORRECTION  TO  THE 

C  ROTATIONAL  ENERGY 

C  DE*J*J*(J+1)**2 

C  WHERE  J   IS  THE  ROTATIONAL  QUAffCUM 

C  NUMBER.      (DEFAULT=0) 

C 

C  BETA THIS  IS  THE  ROTATIONAL  CONSTANT  DUE 

C  TO  THE  CENTRIFUGAL  DISTORTION  AND 

C  FITS  THE  FIRST  TERM  OF  THE  SERIES 

C  BETA*(V+0.5) 

C  (DEFAULT=0) 
C 

C  HV THIS  IS  THE  ROTATIONAL  CONSTANT 

C  ASSOCIATED  WITH  THE  THIRD  ORDER 

C  CORRECTION  TO  THE  ROTATIONAL  ENERGY 

C  HV*(J*(J+1))**3 

C  WHERE  J  IS  THE  ROTATIONAL  QUANTUM 

C  NUMBER.   (DEFAULT=0) 

C 

C  AR THIS  IS  THE  SPIN  COUPLING  CONSTANT 

C  AND  AS  SUCH  IS  IMPORTANT  ONLY  IN 

C  HIGH  RESOLUTION  SPECTRA.  TO  INCLUDE 

C  THE  FUNCTION  ROUTINE  F  MUST  BE 

C  ALTERED.  (DEFAULT=0) 

C 

C   LOWER OOMEGA THE  SAME  AS  OMEGA,  EXCEPT  FOR  THE 

C  LOWER  STATE.  {DEFAULT=0) 
C 

C  AALPHA THE  SAME  AS  ALPHA,  EXCEPT  FOR  THE 

C  LOWER  STATE.  (DEFAULT=0) 
C 

C  BBE THE  SAME  AS  BE,  EXCEPT  FOR  THE 

C  LOWER  STATE.  (DEFAULT=0) 
C 

C  BBETA THE  SAME  AS  BETA,  EXCEPT  FOR  THE 

C  LOWER  STATE.  (DEFAULT=0) 
C 

C  DDE THE  SAME  AS  DE  EXCEPT  FOR  THE 

C  LOWER  STATE.  (DEFAULT=0) 
C 

C   SPEC TE THIS  IS  THE  ENERGY  BETWEEN  THE  MIN- 

C  IMA  IN  THE  POTENTIAL  ENERGY  CURVES 

C  OF  THE  STATES  OF  INTEREST.  (MUST 

C  BE  SUPPLIED) 

C  DNUL THIS  IS  THE  DISSOCIATION  ENERGY 

5  VMIN THE  LOWEST  VIBRATIONAL  STATE  OF 

C  INTEREST.   (DEFAULT=0) 

C  VMAX THE  HIGHEST  VIBRATIONAL  STATE  OF 
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C  INTEREST.  (DEFAULT=0) 

C 

C  VVMIN THE  SAME  AS  VMIN  EXCEPT  FOR  LOWER 

C  STATE.  (DEFAULT=0) 

C 

C  VVMAX THE  SAME  AS  VMAX  EXCEPT  FOR  LOWER 

C  STATE.   (DEFAULT=0) 

C 

C  JMIN LOWEST  ROTATIONAL  STATE  OF  INTEREST. 

C  {DEFAULT=1) 

C 

C  JMAX HIGHEST  ROTATIONAL  STATE  OF  INTEREST. 

C  (DEFAULT=99) 

C 

C   TYPE VERTLG THIS  CONSTANT  DECIDES  THE  TYPE  OF 

C  ROTATIONAL  DISTRIBUTION  TO  BE  USED 

C  VERTLG=0  -  GAUSSIAN  WITH  VAR 

C  GIVING  THE  WIDTH  OF 

C  THE  DISTRIBUTION. 

C  VERTLG>0  -  ARBITRARY  DISTRIBUTION 

C  STORED  IN  PMOD. 

C  VERTLG<0  -  BOLTZMANN  DISTRIBUTION 

C  TEMPERATURE  IS  GIVEN 

C  BY  TEMV. 

C 

C  ROTMAX GIVES  THE  MAXIMUM  ROTATIONAL  VALUE  TO 

C  BE  USED  IN  THE  DISTRIBUTIONS. 

C 

C  VAR AS  DEFINED  ABOVE 

C 

C  JTR VALUE  OF  J  AT  WHICH  TRUNCATION  BEGINS 

C 

C  JCUT VALUE  OF  J  AT  WHICH  DISTRIBUTION  IS 

C  SET  TO  0 

C 

C  TABH THIS  DECIDES  THE  TYPE  OF  VIBRATIONAL 

C  DISTRIBUTION  TO  BE  USED. 

C  TABH  =0   -  SET  OF  ENERGIES  ALONE 

C  DETERMINE  DISTRIBUTION 

C  TABHoO    -  TEMPERATURE  TEMP  (V) 

C  DETERMINES  THE  ENERGIES 

C 

C  TEMP TEMPERATURE  OF  A  GIVEN  VIBRATIONAL  LEVEL 

C 

C  VIN STARTING  VIBRATIONAL  LEVEL  OF  UPPER 

C  STATE 

C 

C  WIN STARTING  VIBRATIONAL  LEVEL  OF  LOWER 

C  STATE 
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C   COLL ECM CENTER-OF-MASS  COLLISION  ENERGY 

C 

C  EXOTH EXOTHERMICITY  OF  THE  REACTION 

C 

C  TEMPV TEMPERATURE  OF  MOLECULE  PRIOR  TO 

C  COLLISION 

C 

C   DISTR FCFFAK FRANCK-CONDON  FACTORS 

C 

C  PMOD SUPPLIED  ROTATIONAL  DISTRIBUTIONS 

C 

C  SWITCH USED  TO  CALCULATE  ROTATIONAL  DIS- 

C  TRIBUTIONS  ONLY  ONE  TIME  PER  PASS 

C 

C  IV PROGRAMMER  SUPPLIED  VIBRATIONAL 

C  DISTRIBUTION 

C 

C  PIV PROGRAMMER  SUPPLIED  VIBRATIONAL 

C  DISTRIBUTION 

C 

C  PSCHW ALTERNATIVE  DISTRIBUTION 

C 

C   PLOT LMIN MINIMUM  WAVELENGTH  OF  INTEREST 

C 

C  LMAX MAXIMUM  WAVELENGTH  OF  INTEREST 

C 

C  KAN NUMBER  OF  STEPS  BETWEEN  LMIN  AND 

C  LMAX 

C 

C  FL AREA  OF  EXPERIMENTAL  SPECTRUM  TO 

C  WHICH  CALCULATION  IS  TO  BE  NORMAL- 

C  IZED 

C 

C  LX THIS  IS  THE  DESIRED  LENGTH  OF  THE 

C  SPECTRA 

C 

C  MSKY THIS  IS  THE  NUMBER  OF  EVENTS/CHANNEL 

C  ANALYZER 

C 

C  AVFL6 THIS  IS  THE  RESOLUTION  OF  SPECTRO- 

C  METER  IN  ANGSTROMS 

C 

C  ITOT THE  NUMBER  OF  SUPPORTING  POINTS  IN 

C  THE  SENSITIVITY  CALCULATION 

C 

C  KL THE  LOWER  CALIBRATION  WAVELENGTH 

C 

C  KH THE  HIGHER  CALIBRATION  WAVELENGTH 
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C  EKL LOWER  WAVELENGTH  LIMIT 

C 

C  EKH HIGHER  WAVELENGTH  LIMIT 

C 

C  ASC LENGTH  OF  SCALE  IN  WAVELENGTHS 

C 

C  DLAMD DISTANCE  OF  SCALE  IN  WAVELENGTHS 

C 

C  ALI WAVELENGTH  OF  CALIBRATION  POINTS 

C 

C  SI CALIBRATION  POINTS 

C 

C  IPLOT DECISION  WHETHER  TO  PLOT  ON  GOULD 

C  PLOTTER 

C  IPLOT=0  -  PLOT  WILL  BE  DRAWN 

Q************************************************************* ******** 
Q********************************************************************* 

DIMENSION  PR(200) ,CTG(1024) ,CT(1024) ,FCFFAK(35,35) , 

#  PM0D(4,100),ALI(100),SI(100).PSCHW(30), 

#  IV(4),PIV(4).T(27),AI(27),AINT(27). 

#  WL(27),0MEGA(12),00MEGA(12),ALPHA(6), 

#  AALPHA(6),TEMP(3),PV(30) 
DOUBLE  PRECISION  OMEGA, OOMEGA.TE 

INTEGER  VMIN,VMAX,VVMIN,VVMAX,VIN,VVIN,V,VV 

REAL  JSKY,LA,LX,LMIN,LMAX 

DATA  PR, FCFFAK, IV, PIV/200*0., 1225*0., 4*0, 4*0./, 

#  PMOD.ALI ,SI,PSCHW/400*0. ,100*1 . ,100*1 . ,30*0./, 

#  PV, OMEGA, 00MEGA/30*0., 12*0., 12*0./, 

#  ALPHA,AALPHA,ITEST,CTU,CTC,CTL/6*0. ,6*0. ,0,0. ,0. ,0./, 

#  BE,BETA,DE,HV,AR/0.,0.,0.,0.,0./, 

#  IPLOT, SWITCH/0.0./ 

NAMELIST  /UPPER/OMEGA,ALPHA,BE,BETA,DE,HV,AR 

#  /LOWER/OOMEGA , AALPHA , BBE , BBETA , DDE 

#  /SPEC/TE,DNULL,VMIN,VMAX,VVMIN,VVMAX,JMIN,JMAX 

#  /TYPE/VERTLG,ROTHMAX,VAR,JTR,JCUT,TABH,TEMP,VIN,VVIN 

#  /COLL/ECM,EXDTH,TEMPV 

#  /DISTR/FCFFAK,PMOD, SWITCH, IV, PIV,PSCHW 

#  /PLOT/LMIN,LMAX,KAN,ITOT,ALI,SI,ALMIN,ALMAX, 

#  MSKY, LX,AUFLG,FL,KH, EKH, ASC, SLAMD, IPLOT 
REA0(5, UPPER) 

READ (5, LOWER) 

READ(5,SPEC) 

READ(5,TYPE) 

READ(5,C0LL) 

READ(5,DISTR) 

READ (5, PLOT) 

CALL  PVIB(PV,VMAX, OMEGA, IV, PIV,TEMPV, SWITCH, PSCHW) 

CHAN=FLOAT(KAN) 
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DEL=(LMAX-LMIN)/CHAN 

VMIN=VMIN+1 

VMAX=VMAX+1 

VVMIN=VVMIN+1 

VVMAX=VVMAX+1 

DO  10  I=VVMIN,VVMAX 

VV=I-1 

GVV=0. 

DO  20  L=1.12 
20  GVV=GVV+00MEGA(L)*(VV+0.5)**L 

BVV=BBE 

DO  30  L=l,6 
30  BVV=BVV+AALPHA(L)*{VV+0.5)**L 

DO  10  I1=VMIN,VMAX 

V=I  1-1 

GV=0 

DO  40  L=l,12 
40  GV=GV+0MEGA(L)*(V+0.5)**L 

EE=(ECH+EX0TH)*8064.5 

IF  (EE.GT.O.)  GOTO  50 

WRITE(5,9000)  V 

GOTO  9999 
50  IF  (TABH.GE.O)  GOTO  60 

TEMV=TEMP(1)*(EE-6V)/EE 

GOTO  70 
60  TEMV=TEMP(V) 
70  TV=TE+GV-GVV 

BV=BE 

DO  80  L=l,6 
80  BV=BV+ALPHA(L)*(V+0.5)**L 

DV=DE+BETA*(V+0.5) 

DVV=DDE+0.5*BBETA*(V+0.5) 

IF  (EE,GT,DNULL)  GOTO  90 

AJMAX=SORT((EE-GV)/BV) 

JJMAX=IFIX(AJMAX) 

GOTO  110 
90  DO  100  J=JMIN,JMAX 

GI=F(J,BV,DV,HV,Y) 

IF  (DNULL+0.5*OMEGA(l)-GV-Gl.LE.O.)  GOTO  100 

JJMAX=J 
100  CONTINUE 
110  IF  (JJMAX.GE.JMAX)  JJMAX=JMAX 

CALL  PROT(PR.ROTMAX,VAR,JJMAX,TEMV,BV,V.VERTLG) 

WRITE(6,9010)  ECM.JJMAX 

DO  10  J=JMIN,JJMAX 

JPL1=J+1 

JMI1=J-1 

Y=AR/BV 
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T(l )=TV+X-FF(JPL1 ,BVV,DVV,HVV) 
T(2)=TV+X-FF{J,BVV,DVV,HVV) 
T(3)=TV+X-FF(JMI1 ,BVV,DVV,HVV) 
DO  115  M=l,3 
115  WL(M)=1.E8/T(M) 

WRITE(6,9020)  (WL(N),N=1 ,3) 

Q********************************************************************* 
Q********************************************************************* 

C  THE  AI(I)'S  WHICH  FOLLOW  ARE  THE  HONL-LONDON  FACTORS  FOR  THE  3-PI-U 
C  TO  3-PI-G  TRANSITIONS.  THESE  MUST  BE  REPLACED  FOR  ANY  OTHER  TYPE 
C  SYMMETRY  TRANSITIONS. 

Q********************************************************************* 

AI(1)=(J+2.)*J/(J+1.) 

AI(2)=(2.*J+1.)/(J*(J+1.)) 

AI(3)=(J+1.)*(J-1.)/J 

FCF=FCFFAK(V+1 ,VV+1 ) 

DO  10  M=l,3 

IF  (T(M).LE.O.)  GOTO  10 

AINT(M)=FCF*T(M)**3*AI(M)*PV(V+1)*PR(J+1)/(J*J+1) 

LA=1.E8/T(M) 

X=(LA-LMIN)/DEL 

IF  (X.LE.CHAN)  GOTO  120 

CTU=CTU+AINT(M) 

GOTO  10 
120  IF  (X.GT.O.)  GOTO  130 

CTL=CTL+AINT(M) 
130  CTC=CTC+AINT(M) 

N=IFIX(X)+1 

CALL  SENS(LA,LMIN,LMAX,S,I  TEST,ALMIN,ALMAX,ITOT,SI) 

CT(N)=CT(N)+AINT(M)*S 
10  CONTINUE 

ANGPK={LMAX-LMIN)/KAN 

BETA=AUFLG/ANGPK 

IBETA=IFIX(BETA) 

IBETA1=IBETA-1 

L=KAN-IBETA 

DO  140  I=IBETA,L 

SUM=1 . 

CTG(I)=CT(I) 

DO  150  J=1,IBETA1 

J1=I-J 

J2=I+J 

CTG(I)=CTG(I)+(IBETA-J)*CT(J1)/IBETA+(IBETA-J)*CT(J2)/IBETA 
150  SUM=SUM+2.*(IBETA-J)/IBETA 
140  CTG(I)=CTG(I)/SUM 

SCTG=0. 

DO  160  N=l.  KAN 
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160  SCTG=SCTG+CT6(N) 

DO  170  N=l,KAiN 
170  CTG(N)=CTG(N)*FL/SCTG 

IF  (IPLOT.NE.O.)  GOTO  9900 

CALL  PL0TS{12.,18.,0,1,1.,2.) 

CALL  AXIS(0.,9.,'WAVELENGTH(A)',-13,12.,0.,ASC,DLAMD) 

DO  180  1=1, KAN 

Y=CTG(I)/MSKY 

X=LX*I/KAN 
180  CALL  PL0T(X,Y,2,2) 

CALL  PL0T(0.,0.,999) 
9900  WRITE  (6, UPPER) 

WRITE(6. LOWER) 

WRITE(6,SPEC) 

WRITE(6,TYPE) 

WRITE(6,TYPE) 

WRITE(6,C0LL) 

WRITE(6,DISTR) 

WRITE(6,PL0T) 
9000  FORMAT  (IX, ' INSUFFICIENT  ENERGY  AT  V  =  ',15) 
9010  FORMAT  (IX.'THE  CENTER-OF-MASS  ENERGY  IS  ',E12.5./. 

#'THE  MAXIMUM  ALLOWED  J  IS  ' ,  15) 
9020  FORMAT  (1X,3E12.5) 
9999  STOP 

END 
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PVIB 

SUBROUTINE  PVIB(PV,VMAX, OMEGA, IV, PIVJEMPV, SWITCH, PSCHW) 

INTEGER  V,VMAX,V1,VMAX1,VSAVE,SWITCH,A 

DIMENSION  PV(1 ) ,0MEGA(1 ) ,IV(1 ) ,PIV(1 )  ,PSCHW(1 ) , 
#        A(50),G(30) 

DATA  IBLANK/'  7,ISTAR/'*7 

VMAX1=VMAX+1 

PPV=0. 

IF  (SWITCH. GT.O)  GOTO  100 

DO  10  V1=1,VMAX1 
10  PV(V1)=PSCHW(V1) 

GOTO  1000 
100  IF  (SWITCH. NE.l)  GOTO  200 

DO  110  V1=1,VMAX1 

G(V1)=0. 

DO  105  L=l,12 
105  G(V1)=G(V1)+0MEGA(L)*(V1-0.5)**L 
110  PV(V1)=EXP(-G1(V1)/TEMPV) 

GOTO  1000 
200  IF  (IV(l).GE.VMAX)  WRITE(6,2000) 

N=IV(1)+1 

DO  290  V1=1,N 
290  PV(V1)=PIV(1) 

DO  300  1=1,3 

IL0=IV(I)+1 

IHI=IV(I+1)+1 

DO  300  V1=IL0.IHI 
300  PV(V1)=(PIV(I+1)-PIV(I))*(VI-(IV(I)+1))/(IV(I+1)-IV(I))+PIV(I) 

IV41=IV(4)+1 

DO  310  V=IV41,VMAX1 

V1=V+1 
310  PV(V1)=0. 
1000  DO  1100  V1=1,VMAX1 
1100  PPV=PPV+PV(V1) 

DO  1200  V1=1,VMAX1 
1200  PV(V1)=PV(V1)/PPV 

VSAVE=1 

DO  1300  V1=1,VMAX1 

IF  (PV(Vl).LT.PV(VSAVE))  GOTO  1300 

VSAVE=V1 
1300  CONTINUE 

DO  1400  V1=1,VMAX1 

V=V1-1 

PL=60*PV(V1)/PV(VSAVE) 

KK=IFIX(PL) 

KK1=KK-1 

DO  1350  K-1,KK1 
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♦  i. 
1350  A(K)=IBLANK 

1400  A(KK)=ISTAR  -  ; 

2000  F0RMAT{1X,'IV(V)  IS  GREATER  THAN  VMAX') 

RETURN 

END 
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PROT 

SUBROUTINE  PROT(PR,ROTMAX,VAR,JTR,JCUT,TEMV,BV,V,BERTLG,PMOO) 

DIMENSION  PR(100),PM0D(3,100) 

INTERGER  V 

REAL  MJ 

SMJ=0. 

SUMG=0. 

N1=JCUT+1 

IF  (VERTLG)  110,10,210 

C     VERTLG=0  YIELDS  GAUSSIAN  DISTRIBUTION  OF  THE  ROTATIONAL  STATES 
C  WITH  VAR  (THE  VARIANCE)  GIVING  THE  WIDTH  OF  THE 

C  DISTRIBUTION 

10  DO  60  L=1,N1 

K=L-1 

IF  (K-JTR)  40.40,50 
40  PR(K+1)=(2*K+1)*EXP(VAR*K) 

GOTO  60 
50  IF  (K=JCUT)  55,55,56 

55  PR(K+1 )=PR(JTR+1 )*(JCUT-K)/(JCUT-JTR) 
GOTO  60 

56  PR(K+1)=0. 

60  SUNG=SUMG+PR(K+1) 
GOTO  1000 

C     VERTLG<0  YIELDS  A  BOLTZMANN  DISTRIBUTION  WITH  THE  TEMPERATURE 
C  OF  THE  DISTRIBUTION  GIVEN  BY  TEMV 

Q************************************ ****************************** 

110  DO  160  L=1,N1 

K=L-1 

IF   (K-JTR)   140.140,150 
140  PR(K+1)=(2*K+1)*EXP(-BV*K*(K+1)/(0.695*TEMV))*BV/TEMV 

GOTO  160 
150  IF   (K-JCUT)   155,155,156 

155  PR(K+1)=PR(JTR+1)*(JCUT  K)/(JCUT-JTR) 
GOTO  160 

156  PR(K+1)=0. 

160  SUMG=SUMG+PR(K+1) 
GOTO  1000 

Q******************************************************************** 

C  VERTLG>0  YIELDS  A  DISTRIBUTION  SET  ARBITRARILY  BY  THE 
C  PROGRAMMER.  THESE  DISTRIBUTIONS  ARE  DEPENDENT  UPON 

C  THE  PARTICULAR  VIBRATIONAL  STATE  BEING  CONSIDERED 

C  AND  ARE  PROVIDED  IN  THE  PMOO  ARRAY 

Q******************************************************************** 

210  DO  220  L=1,N1 
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K=L-1 

PR(K+1 )=(K+1 )*PM0D(V+1 ,K+1 ) 
220  SUMG=SUMG+PR(K+1) 

Q****** ************************************************************** 

C     THE  ROTATIONAL  DISTRIBUTION  IS  NOW  NORMALIZED  TO  A  UNIT  AREA 

Q******************************************************************** 

1000  DO  1010  L=1,N1 

K=L-1 
1010  PR(K+1)=PR(K+1)/SUMG 

RETURN 

END 
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SENS  ,", 

SUBROUTINE  SENS(LA,LMIN,LMAX,S,ITEST,ALMIN,ALMAX,ITOT,ALI ,SI) 
REAL  LA.LMIN.LMAX 

Q******************************************************************** 

C     ALMIN  AND  ALMAX  ARE  THE  WAVELENGTHS  (IN  ANGSTROMS)  FOR  WHICH  \'.^-> 

C  THE  SENSITIVITY  CALIBRATION  IS  AVAILABLE  '":" 

IF  (ALMIN-LMIN)  20,20,30 
20  IF  (ALMAX-LMAX)  40,22,22 
22  ITEST=1 

DO  25  I=1,IT0T 

IF  (LA-ALI(I))  26,26,25 

25  CONTINUE 

26  NX=II-1 
NX1=NX+1 
S=SI(NX)+(SI(NX1)=SI(NX))*(LA-ALI(NX))/(ALI(NX1)-ALI(NX)) 

GOTO  999 
30  PRINT  2100 

GOTO  999 
40  PRINT  2110 
2100  FORMAT  (IX, 'WAVELENGTH  OUT  OF  RANGE  LOW') 
2110  FORMAT  (IX, 'WAVELENGTH  OUT  OF  RANGE  HIGH') 
999  RETURN 
END 


FUNCTION  F(J,BV,DV,HV,Y) 

F=BV*(J*( J+1 )-DV*(J*(J+l ) )**2+HV*( J*( J+1 ) )**3 

RETURN 

END 


FF 

FUNCTION  FF(J,BVV,DVV,HVV) 

FF=BVV*J*( J+1 )-DVV*(J*( J+1 ) )**2+HVV*( J*(J+1 ) )**3 

RETURN 

END 
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